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Abstract 

We study a coupled system of two parabolic equations in one space dimension. This system is singular 
because of the presence of one term with the inverse of the gradient of the solution. Our system describes 
an approximate model of the dynamics of dislocation densities in a bounded channel submitted to an 
exterior applied stress. The system of equations is written on a bounded interval with Dirichlet conditions 
and requires a special attention to the boundary. The proof of existence and uniqueness is done under 
the use of two main tools: a certain comparison principle on the gradient of the solution, and a parabolic 
Kozono-Taniuchi inequality. 
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1 Introduction 

1.1 Setting of the problem 

In this paper, we are concerned in the study of the following singular parabolic system: 



Kt = £K X x H T~p x on I X (0, oo) 

Kx 



[1.1) 



pt = (1 + e)p xx - tk x on /x(0,oo), 

with the initial conditions: 

k(x,0) = k°(x) and p(x, 0) = p°(x), (1.2) 
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and the boundary conditions: 



where 

are fixed constants, and 



k(0,.) =«°(0) and k(1,.) =k°(1), 

p(o,.)=p(i,.) = o, 

e > 0, r£l, 
/:=(0,1) 



(1.3) 



is the open and bounded interval of R. 

The goal is to show the long-time existence and uniqueness of a smooth solution of 
(1.1), (1.2) and (1.3). Our motivation comes from a problem of studying the dynamics 
of dislocation densities in a constrained channel submitted to an exterior applied stress. 
In fact, system (1.1) can be seen as an approximate model of the one described in [12]. 
This approximate model (presented in [12] for e = 0) reads: 



% - £ ®xx 



e+ + e- 

6+ + 9- 



6+ 



on 7 x (0, oo), 
on 7 x (0, oo), 



(1.4) 



with r representing the exterior stress field. System (1.4) can be deduced from (1.1), by 
spatially differentiating (1.1), and by considering 



Px 



9 d 



P , « = P + P , 



(1.5) 



which explains the presence of the factor (1 + e) in the second equation of (1.1). Here 
+ and 6~ represent the densities of the positive and negative dislocations respectively 
(see [25, 16] for a physical study of dislocations). 

The part II of this work will be presented in [18]. There, we will show some kind of 
convergence of the solution (p e , k £ ) as e — > 0. 

1.2 Statement of the main result 

The main result of this paper is: 

Theorem 1.1 (Existence and uniqueness of a solution). Let p°, k° satisfying: 

p , K °eC oo (7"), p°(0)= / > (l) = K (0) = 0, «°(i) = i, (1-6) 

f (1 + e)P°xx = TK °x ° n dl 
\ (1 + £)k° xx = rp\ on dl, 

and 



k x > \p x | on I. 



(1.7) 
(1.8) 
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Then there exists a unique global solution (p, re) of system (1.1), (1.2) and (1.3) satisfying 

(p, k) e C 3+a -^(/ x [0, oo)) n C°°(J x (0, co)), Va€(0,l). (1.9) 
Moreover, this solution also satisfies : 

Kx > \Px\ on / x [0, oo). (1-10) 

Remark 1.2 Conditions (1.7) are natural here. Indeed, the regularity (1.9) of the solu- 
tion of (1.1) with the boundary conditions (1.2) and (1.3) imply in particular (1.7). 

Remark 1.3 Remark that the choice re°(0) = and re°(l) = 1 does not reduce the 
generality of the problem, because equation (1.1) does not see the constants and has the 
following invariance: if (p, re) is a solution, then (Ap, Are) is also a solution for any A G R. 

1.3 Brief review of the literature 

To our knowledge, systems of equations involving the singularity in l/n x as in (1.1) 
has not been directly handled elsewhere in the literature. However, parabolic problems 
involving singular terms have been widely studied in various aspects. Fast diffusion 
equations: 

u t - Au m = 0, < m < 1, 

are examined, for instance, in [5, 7, 8]. These equations are singular at points where 
u = 0. In dimension 1, setting u = v x we get, up to a constant of integration: 

v t - mD™ _1 D n = 

which makes appear a singularity like l/v x . Other class of singular parabolic equations 
are for instance of the form: 

Ut = u xx + -u x , (1.11) 

X 

where b is a certain constant. Such an equation is related to axially symmetric problems 
and also occurs in probability theory (see [6, 23]). An important type of equations that 
can be indirectly related to our system are semilinear parabolic equations: 

u t = Au + |u| p_ y p>l. (1.12) 

Many authors have studied the blow-up phenomena for solutions of the above equation 
(see for instance [24, 13]). Equation (1.12) can be somehow related to the first equation 
of (1.1), but with a singularity of the form 1/re. This can be formally seen if we first 
suppose that u > 0, and then we apply the following change of variables u = 1/v. In 
this case, equation (1.12) becomes: 

2|V^| 2 2 _ v 
v t = Av 1 1 v 2 p , 

v 

and hence if p = 3, we obtain: 

v t = Av -\l + 2\Vv\ 2 ). (1.13) 

Since the solution u of (1.12) may blow-up at a finite time t = T, then v may vanishes at 
t = T, and therefore equation (1.13) faces similar singularity to that of the first equation 
of (1.1), but in terms of the solution v instead of v x . 
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1.4 Strategy of the proof 



The existence and uniqueness is made by using a fixed point argument after a slight 
artificial modification in the denominator k x of the first equation of (1.1) in order to 
avoid dividing by zero. We will first show the short time existence, proving in particular 

that 

Kx (x,t) > Vt 2 W + p£(M) > 7(t) > 0, 

for some well chosen initial data and a suitable function j(t) > 0. The only, but danger- 
ous, inconvenience is that the function 7 depends strongly on \\p X xx(-,t)\\L°°(i), roughly 
speaking: 

7 - -||Pxsx||l°°(/)7- (1-14) 

Let us mention that one of the key points here is that ^ < 1 which somehow linearize 
the first equation of (1.1). Nevertheless, standard Sobolev and Holder estimates for the 
parabolic system (1.1) are not good enough to bound HpxxxIIl 00 !;/) m order to prevent 
7 (and as a consequence k x ) from vanishing. On the contrary, a Sobolev logarithmic 
estimate (see Section 2, the parabolic Kozono-Taniuchi inequality, Theorem 2.13) can be 
used in order to obtain a sharp bound of ||Pxxx||l°°(7) of the form 



|Pxxx||l~(J) < E ^1 + log+ ^ , 



where E is an exponential function in time. This allows, with (1.14), to show that the 
function 7 > does not vanish in finite time. After that, due to some a priori estimates, 
we can prove the global time existence. 



1.5 Organization of the paper 

This paper is organized as follows: in Section 2, we present the tools needed throughout 
this work, this includes a brief recall on the L p , C a and the BMO theory for parabolic 
equations. In Section 3, we show a comparison principle associated to (1.1) that will 
play a crucial rule in the long time existence of the solution as well as the positivity of 
k x . In Section 4, we present a result of short time existence, uniqueness and regularity 
of a solution (p, n) of (1.1). Section 5 is devoted to give some exponential bounds on 
(/>, k). In Section 6, we show a control of the W^' 1 norm of p xxx . In a similar way, we 
show a control of the BMO norm of p xxx in Section 7. In Section 8, we use a parabolic 
Kozono-Taniuchi inequality to control the L°° norm of p xxx . In Section 9, we prove our 
main result: Theorem 1.1. Finally, Sections 10 and 11 are appendices where we present 
the proofs of some technical results. 

2 Tools: theory of parabolic equations 

We start with some basic notations and terminology: 
Abridged notation. 
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• It is the cylinder I x (0,T); / is the closure of /; It is the closure of It] di is the 
boundary of I. 

• ll-llip(n) = ll-llp,n> ^ i s an open set, p > 1. 

• is the lateral boundary of It, or more precisely, St = di x (0, T). 

• <9 p 7r is the parabolic boundary of It, i.e. 8 p It = St U (/ x {i = 0}). 

• -D^u = it is a function depending on the parameter y, s£ff. 

• [/] is the floor part of I e R. 

• Q r = Qr(xo,to) is the lower parabolic cylinder given by: 

Q r = {(x,t); \x - x \ <r, t - r 2 < t < t }, r > 0, (x ,t ) G ir- 

• |0| is the n-dimensional Lebesgue measure of the open set SI C R n . 

• ms)(«) = f n u is the average integral of the u G over !]cl n . 
2.1 L p and C Q theory of parabolic equations 

A major part of this work deals with the following typical problem in parabolic theory: 

u t = eu xx + / on I T 

u(x,0) = <f> on / (2.1) 

« = on dIx(0,T), 

where T > and e > 0. A wide literature on the existence and uniqueness of solutions 
of (2.1) in different function spaces could be found for instance in [21], [11] and [22]. We 
will deal mainly with two types of spaces: 

The Sobolev space W^^t), 1 < p < oo which is the Banach space consisting of the 
elements in L p (It) having generalized derivatives of the form D r t D s x u, with r and s two 
non-negative integers satisfying the inequality 2r + s < 2, also in L p (It)- The norm in 
this space is defined as \\u\\ w 2,i {It) = ELoE 2 r+s=i \\ D t DS x u \\p,i T - 

The Holder spaces C e (I) and C^/ 2 (It), £ > a nonintegral positive number. We do 
not recall the definition of the space C (I) which is very standard. The Holder space 
C^/ 2 (7t) is the Banach space of functions v(x,t) that are continuous in It, together 
with all derivatives of the form D r t D x v for 2r + s < £, and have a finite norm \v\j^ = 
Wg+EiloMg, where 

2r+s=j 
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and 



with 



2 

T 

2r+s=[t] 0<£-2r-s<2 



(v)^ = mf{c; \v(x,t) -v(x',t)\ < c\x-x'\ a , (x, t), (x' , t) G I T }, < a < 1, 

(v)^ = inf{c; |u(x,t) -v(x,t')l < c|* — (x,t),{x,t') eh}, < a < 1. 

The above definitions could be found in details in [21, Section 1]. Now, we write down 
the compatibility conditions of order and 1. These compatibility conditions concern 
the given data 4>, <E> and / of problem (2.1). 

Compatibility condition of order 0. Let <fi G C{I) and <1> G C(St)- We say that the 
compatibility condition of order is satisfied if 

*\ai = $ Lo- (2-2) 

Compatibility condition of order 1. Let (f> G C 2 {I), <I> G C 1 ^) and / G C(7r). 
We say that the compatibility condition of order 1 is satisfied if (2.2) is satisfied and in 
addition we have: 

{e^ x + f)\ 9I = — t=Q (2.3) 

We state two results of existence and uniqueness adapted to our special problem. We 
begin by presenting the solvability of parabolic equations in Holder spaces. 

Theorem 2.1 (Solvability in Holder spaces, [21, Theorem 5.2]). Suppose 
< a < 2, a non-integral number. Then for any 

</> G C 2+a {I), $ G C 1+a/2 (Sr~) and f G C a ' a/2 (h) 

satisfying the compatibility condition of order 1 (see (2.2) and (2.3)), problem (2.1) has 
a unique solution u G C 2+a ' 1+a / 2 (/x) satisfying the following inequality: 

Ht a) < c H (i/ig> + \<t>\f +a) + mt a/2) ) > ( 2 - 4 ) 

for some c H = c H (e,a,T) > 0. 

Remark 2.2 (Estimating the term c H of (2-4))- The constant appearing in the 
above Holder estimate (2.4) can be estimated, using some time iteration, as c H {e,a,T) < 
e c ( T+1 \ where c = c(e,a) > is a positive constant. 

We now present the solvability in Sobolev spaces. Recall the norm of fractional Sobolev 
spaces. If / G W£(a, b), s > and 1 < p < oo, then 

W-{a,b) ~ WfW W l%, b) +\J a J a \ x _ y \l+(s-[s]) P I • ( 2 - 5 ) 
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Theorem 2.3 (Solvability in Sobolev spaces, [21, Theorem 9.1]). Let p > 1, 

e > and T > 0. For any f G LP(I T ), (f> G Wp' 2/p (I) and $ G w£~ 1/?p (,St), wtt 
p / 3/2 (p = 3/2 is called the singular index) satisfying in the case p > 3/2 the com- 
patibility condition of order zero (see (2.2)), there exists a unique solution u G W 2,1 (It) 
of (2.1) satisfying the following estimate: 

W U Wwt>\l T ) ^ C (WfWpJT + U\\ w 2-2/ P(I) + \\H W 1-V*P {St) ) , (2-6) 

for some c = c(e,p,T) > 0. 

For a better understanding of the spaces stated in the above two theorems, especially 
fractional Sobolev spaces, we send the reader to [1] or [21]. The dependence of the 
constant c of Theorem 2.3 on the variable T will be of notable importance and this what 
is emphasized by the next lemma. 

Lemma 2.4 (The constant c given by (2.6): case = and = 0). Under the 
same hypothesis of Theorem 2.3, with = and = 0, the estimate (2.6) can be written 
as: 



H 7=p" + \\uxx\\ P ,i T + IKIU/t ^ c||/|| Pi 7 T , (2.7) 



T y/f 

where c = c(e,p) > is a positive constant depending only on p and e. 

The proof of this lemma will be done in Appendix A. Moreover, We will frequently make 
use of the following two lemmas also depicted from [21]. 

Lemma 2.5 (Sobolev embedding in Holder spaces, [21, Lemma 3.3]). 

(i) (Case p > 3). For any function u G Wp' 1 (It), if a = I — 3/p > 0, i.e. p > 3, then 
u G C 1+a '~^~ (It) with |u|^ +a ^ < c||u|| w 2,i^.^, c = c(p,T) > 0. However, in terms of 
u x , we have that u x G C a ' Q / 2 (/r) satisfies the following estimates: 

IM|oo,/ T <c{6 a (\\ut\\p,i T + \\u xx \\pj T ) + 5 a ~ 2 \\u\\ p j T }, c = c(p) > 0, (2.8) 

(ui))[y < c|||u t || Pi / T + Hu^IIp./t + ^ IMk/ T j . c = c (p) > 0. 

(ii) (Case p > 3/2). If u G Wp^^x) with p > 3/2, then u G C(It), and we have the 
following estimate: 

IM|oo,/ T < c {<5 2 - 3/p (|M| Pi / T + ||«xx|| P) / T ) + <5" 3/p IMIp,/t} . c = c(p) > 0. (2.9) 
In the above two cases 5 = min{l/2, VT}. 



Lemma 2.6 (Trace of functions in W 2,1 (It), [21, Lemma 3-4])- Ifu€ W 2,1 (It), 
p>l, then for2r + s<2- 2/p, we have D r t D s x u\ t=() G Wp~ 2r ~ s ~ 2/p {I) with 

||«|| M ,2- 2r -.- 2/P(J) < c{T)\\u\\ w 2,i {It) . 
In addition, for 2r + s < 2 — l/p, we have DlD x u\-g^ G Wp r s ^ 2 1 ^ 2p (St) with 



|«|| Wp l-r-./2-l/2p (Sf ) < c(T)||«|| w a,l (jT) . 
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A useful technical lemma will now be presented. The proof of this lemma will be done 
in Appendix A. 



Lemma 2.7 (L°° control of the spatial derivative). Let p > 3 and let < T < 1/4 

(this condition is taken for simplification). Then for every u £ Wp'^ilx) with u = on 
d p (Ir) in the trace sense (see Lemma (2.6)), there exists a constant c(T,p) > such that 

IKHoo,/ T < c(T,p)\\u\\ w 2,iy T y with c(T,p) = c(p)T 2 p -> as T -> 0. 
2.2 BMO theory for parabolic equation 

A very useful tool in this paper is the limit case of the LP theory, 1 < p < oo, for parabolic 
equations, which is the BMO theory. Roughly speaking, if the function / appearing in 
(2.1) is in the LP space for some 1 < p < oo, then we expect our solution u to have 
ut and u xx also in L v . This is no longer valid in the limit case, i.e. when p = oo. In 
this case, it is shown that the solution u of the parabolic equation have Ut and u xx in 
the parabolic/anisotropic BMO (bounded mean oscillation) space that is convenient to 
present some of its related theories. 

Definition 2.8 (Parabolic/ Anisotropic BMO spaces). A function u € L] oc (It) is 
said to be of bounded mean oscillation, u £ BMO(Lt), if the quantity 



sup ( 1777 / \ u - m Qr( u )\) 

2rC/ T \ \Qr\ jQ r / 



is finite. Here the supremum is taken over all parabolic lower cylinders Q r . 

Remark 2.9 The parabolic BMO(Lt) space, which will be refereed, for simplicity, as 
the BMO(Lt) space, and sometimes, where there is no confusion, as BMO space, is a 
Banach space (whose elements are defined up to an additive constant) equipped with the 
norm 

\\u\\bmo(i t ) = sup (7777/ \u-m Qr (u)\). 

We move now to the two main theorems of this subsection, the BMO theory for parabolic 
equations, and the Kozono-Taniuchi parabolic type inequality. To be more precise, we 
have the following: 

Theorem 2.10 (BMO theory for parabolic equations in the periodic case). 

Take < T\ < T. Consider the following Cauchy problem: 

(u t = eu xx + f on Rx(0,T), 
[ u(x,0) =0. 

If f £ L°°(1R x (0, T)) and f is a 2L-periodic function in space, i.e. f{x + 2,t) = 
f(x,t), then there exists a unique solution u G BMO(M. x (0,T)) of (2.10) with 
Ut, u xx S BMO(M. x (0, T)). Moreover, there exists c > that may depend on T\ 
but independent of T such that: 

IKIIbMO(1Rx(0,T)) + \\ u xx\\bMOQZLx(0,T)) < c [||/IIbMO(IRx(0,T)) + m 2Ix{Q,T) (l/l)] ■ (2- 11 ) 
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The proof of this theorem will be presented in Appendix B. Our next tool (see Theorem 
2.13) shows an estimate involving parabolic BMO spaces. This estimate is a control of 
the L°° norm of a given function by its BMO norm and the logarithm of its norm in a 
certain Sobolev space. It can also be considered as the parabolic version on a bounded 
domain It of the Kozono-Taniuchi inequality (see [20]) that we recall here. 

Theorem 2.11 (The Kozono-Taniuchi inequality in the elliptic case, [20, The- 
orem 1]). Let 1 < p < oo and let s > n/p. There is a constant C = C(n,p,s) such 
that, for all f G Wp(M. n ), the following estimate holds: 

II/IU.R" < C (l + ||/|| B MOe(R») (l+log + ||/||^(R"))) , log+ = max(0,log). (2.12) 

Remark 2.12 It is worth mentioning that the BMO e norm appearing in (2.12) is the 
elliptic BMO e norm, i.e. the one where the supremum is taken over ordinary balls 

B r {X ) = {X G R n ; \X - X \ < r}. 

The original type of the logarithmic Sobolev inequality was found in [3, 4] (see also [9]), 
where the authors investigated the relation between L°°, and W£ and proved that 
there holds the embedding 

(l + log"- 1 (l + IKHw^CM-))) , sp>n 

provided IHI^/fc < 1 for kr = n. This estimate was applied to prove existence of global 
solutions to the nonlinear Schrodinger equation (see [3, 14]). 

In our work, we need to have an estimate similar to (2.12), but for the parabolic 
BMO space and on the bounded domain It- This will be essential, on one hand, to 
show a suitable positive lower bound of k x (k given by Theorem 1.1), and on the other 
hand, to show the long time existence of our solution. Indeed, there is a similar inequality 
and this is what will be illustrated by the next theorem. 

Theorem 2.13 (A parabolic Kozono-Taniuchi inequality, [17, Appendix B2], 
[19]). Let v G W 2 ' (It), then there exists a constant c = c(T) > such that the estimate 
holds 

\\v\\oo,i T < c\Hbmo(i t ) i 1 + log+ IHWVt)) ' ( 2 - 13 ) 
where BMO(I T ) = BMO{I T ) n L X {I T ), and for v £ BMO(I T ), 

W v \\bmo(i t ) = \\ v \\bmo(i t ) + IMIli(/ t )- 
This inequality is first shown over R x xR t , then it is deduced over It- 

3 A comparison principle 

Proposition 3.1 (A comparison principle for system (1.1)). Let 

(p, k) G (C 3+a '^ (ir) f for some < a < 1, 
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be a solution of (1.1), (1-2) and (1.3) with k x > 0, and the initial conditions p° , k° 
satisfying: 

4 >\/lo + (P°x? on I, 7 o G (0,1). (3.1) 
Choose (5 = (5{e, r) > large enough. Let the function j(t) satisfies: 

< - (CO + \\Pxxx(;t)\\ L °°(I)) , co = co(£,(3,t), 
7 (0) = 70 /2. 



De/me M(a;,t) := cosh(f3x){n x (x, t) - yVW + {p x (x,t)) 2 } for (x,t) G J T and T > 0. 
T/ien m(i) := min M(x,i) satisfies m(t) > j 2 (t) for all t G [0, T]. In particular, we have 

x£l 



K x (x,t) > vV(*) + p x (x,t), in I T . (3.3) 
Proof. Throughout the proof, we will extensively use the following notation: 



G a {y) = Va 2 + y 2 , a,yeR. 

Without loss of generality (up to a change of variables in (x, t) and a re-definition of r), 
assume in the proof that 

J = (-1,1). 

Define the quantity M by: 

M(x,t) = K x (x,t) - G y( t)(p x (x,t)), (x,t) G I T , 
7(t) > is a function to be determined. The proof could be divided into five steps. 
Step 1. (Partial differential inequality satisfied by M) 
We first do the following computations on It'- 

M t = n xt - d^ Px ) Pxt - J2_ (3.4) 

V 7 + Px 

M x = k xx - G y (p x )p xx , M xx = k xxx - G^(p x )p 2 xx - G^(p x )p xxx . (3.5) 
Deriving (1.1) with respect to x, we deduce that 

2 

_ , Pxx | PxPxxx PxPxx^xx 

l^xt — £K-xxx H I 2 T~p xx , 

K x K x K x (,3.DJ 

Pxt = (1 + &)Pxxx ~ TK XX . 



We set 



r = ^===, F y (y) = y -7arctan(y/7). 
V 7 + P£ 
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Doing again some direct computations, and using (3.4), (3.5) and (3.6), we obtain 
M t = eM xx + (tG' j{Px ) - ^ W + - M 



K x / \ K x 



„2 

+ EG"{ Px )pl x + !2f[Gi(p x ) - G'j Px ) Px ] - r(l - F{ Px )) Pxx - V. 



(3.7) 



K x 



Using Young's inequality 2ab < a 2 + b 2 , we have: 

TJ 2 \ P xx\ £J 2 plx l 2 ^ 2 ( -v 

7 2 + Pi ~ (7 2 + ^) 3/2 4 ev /^T73T 1 • j 

Plugging (3.8) into (3.7), and using some properties of G 1 and F 1 , we get: 

M t > eM xx + (tG' 7 ( Px ) - E^) M x+ (*f- M 



77 



4e V / 7 2 + P 2 vV + Pi 
Step 2. (The boundary conditions for M) 

The boundary conditions (1.3), and the PDEs of system (1.1) imply the following equal- 
ities on the boundary (using the smoothness of the solution up to the boundary), 

£K X x + -t Px = on dl x [0, T] 

K x (3.9) 

^l+e)p xx -TK x = on 8Ix[0,T]. 

In particular (3.9) implies 

M x = —^-G\( Px )M on dIx[0,T}. (3.10) 

To deal with the boundary condition (3.10), we now introduce the following change of 
unknown function: 

M(x, t) = cosh((3x)M(x, t), (x, t) £ It- 
We calculate M on the boundary of I to get: 

M x =(ptaj&(J3x)-j^Gf^(p x )jM on dl x [0,T]. (3.11) 

We claim that, for any fixed time t, it is impossible for M to have a positive minimum 
at the boundary of /. Indeed we have the following two cases: 

M has a positive minimum at x = 1 =>■ M x < 0; 

M has a positive minimum at x = — 1 =>• M x > 0. 
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Both cases violate the equation (3.11) in the case of the choice of (3 = (3(e, r) large enough, 
and hence the minimum of M is attained inside the interval /. Direct computations give: 



M t > eM xx + 



tG'J Px ) - - 2/fetanh(/3x) 



K x 



— cosh(/?x)7 2 r 2 cosh(/3x)77 
• ■ •'• . ■-— - — — j 

X 



4e^ 7 2 + pi VT+ff 



+ 



p 2 xx PxxxG 7 (p x ) 



-/?tanh(/3x) tGJ Px ) - + ef3 z {2 tanh 2 (/?x) - 1) 



X 



M. 
(3.12) 



Step 3. (The inequality satisfied by the minimum of M) 

Let 

m(t) = min M(x,i). 

X&I 

Since the minimum is attained inside /, and since M is regular, there exists xo{t) £ / 
such that m(t) = M(xo(t),t). We remark that we have: 

M x (x (t),t) = 0, and M xx (x (t) , t) > 0, 

and hence, using (3.12), we can write down the equation satisfied by m, we get (indeed 
in the viscosity sense): 

R 

/ ; ^ 

^ ( p xx PxxxGy(p x ) 

71 " > J ^—vm-v \ ' ^J\HXJ ~2 

\ ,v x K x \ K x 

cosh(/3x)7 2 r 2 cosh((3x)j^' _ . . 

4e V7 2 + Pi V7 2 + Pi 

(3.13) 

Step 4. (Estimate of the term R) 

We turn our attention now to the term R from (3.13). Using elementary identities, we 
get 



_ > ,Hxx _ i^x, _ /3tanh(/3x) tG ( px ) _ ™i + e/ 3 2 (2tanh 2 (/3x) - 1) m 



R > - P - G '^ - ^pM P l _ ^ {Gj{px)r _ (3.14) 
By (3.1), we know that 

m(0) > 7 2 (0), 

and the continuity of m preserves its positivity at least for short time. Then, as long as 
m is positive, we have 

«* > vV + 7 X - (3-15) 

Let 

= \\Pxxx(;t)\\oo,I- 

By using (3.15) and some basic identities, inequality (3.14) implies: 

R>~ r4-^ ~ Cl > Cl = T + £ + £/?2 - (3 - 16) 
VT + Pl 4 8e 
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Step 5. (The choice of 7 and conclusion) 

When 7' < 0, we deduce from (3.13) and (3.16) that the function m is a viscosity 
super-solution of: 



m t 



+ c\ m 



C27 



77 



\Ar 2 + pI vV + pI ' 



C2 



r 2 cosh /? 
4e 



(3.17) 



We remind the reader that p x = p x (xo(t),t). Take the function 7 satisfying: 



( 1 
7_ 

7 



<-(c + c), c = min(ci,c 2 ), 



7(0) = 7o/2 



Plug m = 7 2 into (3.17), we directly deduce that 7 is a viscosity sub-solution of (3.17), 
and the result follows by comparison. □ 



4 Short time existence, uniqueness, and regularity 

In this section, we will prove a result of short time existence, uniqueness and regularity 
of a solution of problem (1.1), (1.2) and (1.3). 



4.1 Short-time existence and uniqueness of a truncated system 

We denote 

Ia,b '■= I x (0,0 + b), a,b > 0. 
Fix To > 0. Consider the following system defined on It ,t by: 

Kt = £K XX H Tp x on I TotT 

k x 

p t = (1 + e)p xx - tk x on I To:T , 

with the initial conditions: 

k(x,T ) = k Tq (x) and p(x, T ) = p T °(x), (4.2) 
and the boundary conditions: 



k(0, .) = and k(1, .) = 1 for T < t < T + T 
p(0, .) = p(l, .) = 0, for T < t < T + T. 



(4.1) 



(4.3) 



Remark 4.1 (The terms p and a). In all what follows, and unless otherwise precised, 
the terms p and a £ (0, 1) are two fixed positive real numbers such that 

p > 3 and a = 1 — 3/p. 

Concerning system (4.1), (4.2) and (4.3), we have the following existence and uniqueness 
result. 
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Proposition 4.2 (Short time existence and uniqueness). Let p > 3, and To > 0. 

Let 

p To ,K To € C°°(/x {T }) 

be two given functions such that p T °(0) = p T °(l) = k t °(0) = 0, and k t °(1) = 1. Suppose 
furthermore that 

4° > 7o on I x {t = T }, 

and 

WiD'y^D^Wooj <M on Ix{t = T }, a = 1,2, 
where 70 > and Mq > are two given positive real numbers. Then there exists 

T = T* = T* (M , 70, e,T,p) > 0, (4.4) 

such that the system (4-1), (4-%) an d (4-3) admits a unique solution 

(p, K )e(W^(l To , T )f. 

Moreover, this solution satisfies 



k x > 7o/2 on I TotT , (4.5) 

and 

\p x \ < 2M on 7t~^. (4.6) 

Proof. The short time existence is done by using a fixed point argument. Since we are 
looking for solutions satisfying (4.5) and (4.6), we artificially modify (4.1), and look for 
a solution of 

Kt = £K xx + - JTyTT-Tpx m I To ,T . A 

(70/2) + (k x - 7o/2)+ (4.7) 
. Pt = (1 + e)Pxx ~ tk x in I To>T , 

with the truncation function Tq{x) = ,£)+Cl{x>c} — Cl{x<-C}> C > 0, and satisfying 
the same initial and boundary data (4.2), (4.3). Denote 

y = w^(i To , T ). 

For any constant A > 0, let us define D p x and D x as the two closed subsets of Y given 
by: 

D{ = {ue Y; \\u x \\ p j TotT < A, u = p T ° on d p I To ,T} 

and 

T>£ = {« € Y; \\v x \\ pjTo T < A, v = k To on c> p /t ,t}. 

We choose A large enough such that these sets are nonempty. Define the application \£ 
by: 

:D p x x D^i — ► D{x Dl 

{p,k)\ — ► *(p,K) = (p,/c), 
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where (p, k) is a solution of the following system: 

K t = EK xx + - -— ~Tp x m I TT , 

(70/2) + - 7o/2)+ (4.8) 
, ft = (1 + e)Pxx - tk x in J TojT , 

with the same initial and boundary conditions given by (4.2) and (4.3) respectively. 
The existence of the solution of (4.8), (4.2) and (4.3) is a direct consequence of Theo- 
rem 2.3. Taking p(x,t) = p(x,t) — p T °(x) and R(x,t) = n(x,t) — k t °(x), we can easily 
check that (p, R) satisfies a parabolic system similar to (4.8) with (p, R) = on 8 p It ,t- 
Using Sobolev estimates for parabolic equations to the system satisfied by (p,R), par- 
ticularly (2.7), we deduce that for sufficiently small T > 0, we have ||/9 x ||p,/ To t < A, 
||^a;||p,/T t — A, and hence the application ^ is well defined. 

The application $ is a contraction map. Let ^(p, k) = (p, k) and ^>(p', k') = (//, k!). 
Direct computations, using in particular (2.7), give: 

\\p-p'\\y < cVr\\k- k'\\ Y , (4.9) 

and 

\\k - k'\\ y < c\\F\\ p j TqT , (4.10) 

with the function F satisfying: 

Ai A 2 

■ : — " > ^~ r — * 

j?_ s _r(z_i{\ _ T 2Mo {px) , / s , Pxx( T 2M (px) ~ T 2 M (p' x )) 

+ {P 9 )x ~ (7o/2) + (k x - 70 /2)+ {Pxx P ™ )+ (70/2) + (k x - 70 /2)+ 

A 3 



+ Pxx T 2M {p'x) 



1 1 



. (70/2) + (k x - 70 /2)+ (70/2) + (k' x - 70/2)- 

(4.11) 

In order to prove the contraction for some small T > 0, we need to estimate all the terms 
appearing in (4.11). The term A\ can be easily handled. However, for the term A 2 , we 
proceed as follows. We apply the L°° control of the spatial derivative (see Lemma 2.7) 
to the function p — p', we get: 

II (P ~ P%\\oc,i To , T <cT-W \\p- p'\\ Y . (4.12) 
For the term p xx , we apply (2.7), and hence we deduce that 

\\p' xx \\ p ,i To , T <c(M + \). (4.13) 
From (4.12) and (4.13), we deduce that 

7o 

The term A3 could be treated in a similar way as the term A 2 . The above ar- 
guments, particularly (4.9) and (4.10), give the contraction of ^ for small time 
T = T*(Mo,7o,£,T,p) > 0. Finally, inequalities (4.5) and (4.6) directly follow using 
the Sobolev embedding in Holder spaces (Lemma 2.5). □ 
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4.2 Regularity of the solution 

This subsection is devoted to show that the solution of (4.1), (4.2) and (4.3) enjoys 
more regularity than the one indicated in Proposition 4.2. This will be done using a 
special bootstrap argument, together with the Holder regularity of solutions of parabolic 
equations. 

Proposition 4.3 (Regularity of the solution: bootstrap argument). Under the 
same hypothesis of Proposition 4-2, let p T ° and k t ° satisfy: 

{ (1 + e)pl\ = tkI° at dl, 

{ y r t (4.14) 

\ (1 + e)«£°= Tp? at dl. 

Then the unique solution (p, k) given by Proposition 4-2 is in fact more regular. Precisely, 
it satisfies for a = 1 — 3/p: 

p, K e C 3+a ^(h~^) n C°°(I x (T , T + T)), (4.15) 

where T is the time given by Proposition 4-2. 

Proof. For the sake of simplicity, let us suppose that To = 0. 

The Holder regularity. Since k 6 Wp^fa), we use Lemma 2.5 to deduce that k x 6 
C a > a / 2 (I T ). We apply the Holder theory for parabolic equations Theorem 2.1, to the 
second equation of (4.1) (using in particular the regularity of the initial data p°), we 
deduce that: 

P eC 2+a ' 1+a/2 (h). (4.16) 

Here the compatibility condition is satisfied by (4.14). Using (4.16) and (4.5), we deduce 
that px P xx — rp x £ C a ' a ^ 2 (lT) and similar arguments as above give that: 

K eC 2+a ' 1+a/2 (h). (4.17) 

Repeating the above arguments, using this time (see (4.17)) that k x G C 1+a '~^~ (It), 
and hence 

p( z C 3+a ^{h), (4.18) 
where (4.18) directly implies that px £ xx — rp x G C 1+a '^~ (It), and therefore 

Ke C 3+a '^(h). (4.19) 

The compatibility condition of order 1 which is needed to apply Theorem 2.1 is always 
satisfied by (4.14). The Holder regularity of (p, k) directly follows from (4.18) and (4.19). 

The C°° regularity. In order to get the C°° regularity, we argue as in the case of 
the Holder regularity (bootstrap argument). In this case the compatibility condition is 
replaced by multiplying by a test function that vanishes near t = 0. □. 
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5 Exponential bounds 



In this section, we will give some exponential bounds of the solution given by Proposition 
4.2, and having the regularity shown by Proposition 4.3. It is very important, throughout 
all this section, to precise our notation concerning the constants that may certainly vary 
from line to line. Let us mention that a constant depending on time will be denoted by 
c(T). Those which do not depend on T will be simply denoted by c. In all other cases, 
we will follow the changing of the constants in a precise manner. 

Proposition 5.1 (Exponential bound in time for p x and n x ). Let 

p,K& C 3+a >'^r{I x [0, oo)) n C°°(J x (0, oo)), 

be a solution of (1.1), (1.2) and (1.3), with p°(0) = p°(l) = 0, k°(0) = and k°(1) = 1. 
Suppose furthermore that the function 

Px 

B = — satisfies ||5||l~(/x(o,oo)) < L 

Kx 

Then, for small T* = T*(e,T,p) > 0, and A = 1 + ||p || VF 2-2/ P ^ + ||«°|| w -2-2/ P ^, we 
have for all t > 0: 

\ Px \^,\n x \^<cAe ct , (5.1) 
and c is a fixed constant independent of the initial data. 

Proof. We use the special coupling of the system (1.1) to find our a priori estimate. 
Roughly speaking, the fact that k x appears as a source term in the second equation of 
system (1.1) permits, by the L p theory for parabolic equations, to have L p bounds, in 
terms of ||k x |Ip,/t5 on Ar an< ^ Pxx which in their turn appear in the source terms of the 
first equation of (1.1) satisfied by k. All this permit to deduce our estimates. To be more 
precise, let T > an arbitrarily fixed time, the proof is divided into four steps: 

Step 1. (estimating k x in the L p norm) 

Let k be the solution of the following equation: 

on It 

(5.2) 

on d p I T - 

As a solution of a parabolic equation, we use the L p parabolic estimate (2.6) to the 
function k to deduce that: 

ll«'ll^(7 T )<^)(ll«°ll w -/ P(J) + l), (5-3) 

where the term 1 comes from the value of n' = k on St- Take 

K = K — K , (5-4) 
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then the system satisfied by k reads: 

, i f . Px Pxx 7- 

K t = K xx - (K t - SK XX ) H Tp x Oil 1 T 

n x 

R = on d p I T . 

Using the special version (2.7) of the parabolic L p estimate to the function k, we obtain: 

II^2:IIp,/t — cv^r ^||^tllp,/ T + II^xxIIp^t Wpxx \\p,ir \\p% Wpjr^j ' (5-5) 

where we have plugged into the constant c the terms e, r, p and ||-B||oo- Combining (5.3), 
(5.4) and (5.5), we get: 

\K\\ P ,i T < c(T) (||« || w 2-2/p (J) + l) + cVT\\p\\ w 2,i {lT) . (5.6) 

The term ||p|| w 2,i( jT ) appearing in the previous inequality is going to be estimated in 
the next step. 

2 1 

Step 2. (estimating p in the W p ' norm) 

As in Step 1, let p , p be the two functions defined similarly as k , R respectively (see 
(5.2) and (5.4)). The function p satisfies an inequality similar to (5.3) that reads: 

\\p'Ww1>\i t ) ^ c ( T )\\p°\\ w *- 2/p (iy ( 5 - 7 ) 

The term 1 disappeared here because p = p = on St- We write the system satisfied 
by p, we obtain: 

ipt = {l + e)p xx + ((1 +e)p xx - p t )- tk x on I T 
\p(x,0)=0 on d p I T , 

hence the following estimate on p, due to the special L p interior estimate (2.7), holds: 

\\p\\ w ^{i T ) - c (ll^tllp^T + Wp'xxWpJt + IMIp,/t) • ( 5 - 8 ) 

Again, we have plugged e, r and p into the constant c, and we have assumed that T < 1. 
Combining (5.7) and (5.8), we get in terms of p: 

<c(T)\\ P \\ w ^- 2 /p( I ) + c \\ K x\\p,i T - (5.9) 

We will use this estimate in order to have a control on ||« x || P) j T for sufficiently small time. 

Step 3. (Estimate on a small time interval) 

From (5.6) and (5.9), we deduce that: 

||«x||p,J T < c(T) (ll«°ll w -2-2/ P(7) + ||p ||^2-2/ P(/) + l) + cVt\\k x \\ pJt . (5.10) 

18 



Let us remind the reader that all constants c and c(T) have been changing from line to 
line. In fact, the important thing is whether they depend on T or not. Let 

T* = ^-3, c is the constant appearing in (5.10), 

we deduce, from (5.10), that 

||«x||p,/ x * < C3 ^||k || h/ ^-2/ P( . / j + ||p°|| w .2-2/ P ^ + 1^ , 

where C3 = C3(T*) > is a positive constant which depends on T* . Recall the special 
coupling of system (1.1), together with the above estimate, we can deduce that: 

ll(P> K )Wwl>\i T *) ^ c 4 (\\ K °\\ W ?- 2/p (i) + \\P°\\w*- 2/p (i) + X ) ' ( 5 - U ) 

with C4 = C4(T*) > is also a positive constant depending on T* but independent of the 
initial data. 

Step 4. (The exponential estimate by iteration) 

Now we move to show the exponential bound. Set 

f{t) = ||(p,K)|| w -a.i (Jx(t>t+T . ))) and g(t) = \\k(; t)\\ w 2-a/ P (/) + \\p(-,t)\\ w 2-2/ P{iy 

Using estimate (5.11) of Lemma 2.6, together with estimate (5.11) of Step 3, we get 

g{T*) < c 5 /(0) < c 5 c 4 (g(0) + 1), c 5 = c 5 (T*). 

In this case, the Sobolev embedding in Holder spaces (see Lemma 2.6), and the time 
iteration give immediately the result. □ 

Proposition 5.2 (Exponential bound in time for p xx ). Under the same hypothesis 
of Proposition 5.1, and for some T* = T*(e,T,p) > 0, we have: 

|p|g+? < cAe c \ t>0, (5.12) 



1 ll 00, 4- II AT II n - • _L I o0|( 2 +") 

J T ,,2-2/p /rN -+- \\K L.,2 

constant independent of the initial data 



where A = 1 + ||p || H/ 2-2/ P ^^ + ||«; || w 2-2/p^ + |p°|j , and c > is a fixed positive 



Proof. The proof is very similar to the proof of Proposition 5.1. It uses in particular 
the Holder estimate for parabolic equations (namely (2.4)), the Holder embedding in 
Sobolev spaces (Lemma 2.5), and finally the iteration in time. □ 

Remark 5.3 We can not obtain, using similar arguments as in the proof of Proposition 
5.2, a similar exponential bound (5.12) for the term H^^- This is due to the presence 
of the term 1/k x in the equation involving k. 
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2 1 

6 An upper bound for the W 2 ' norm of p. 



2 1 

This section is devoted to give a suitable upper bound for the W 2 ' norm of p xxx . This 

2 1 

result will be a consequence of the control of the W 2 ' norm of K t and k xx . The goal 
is to use this upper bound in the parabolic Kozono-Taniuchi inequality (see inequality 
(2.13) of Theorem 2.13) in order to control the L°° norm of p xxx . 

Let us fix T\ > 0. In this section (Section 6) and in the following section (Section 7), 
we will obtain some estimates on the solution (p, n) on the time interval (0, T), with 

T > Ti > 0. (6.1) 

In these estimates, we will precise the dependence on T which involves some constants 
depending on T\ that may blow up as T\ goes to zero. Consider the following hypothesis: 

(HI) The function k x satisfies: 

n x {x,t) > 7(t) > 0, 
where 7 (i) is a positive decreasing function with 7(0) 

Let 

V = I T , 

we start with the following lemma. 

Lemma 6.1 (W^ 1 bound for n t and k xx ). Under hypothesis (HI), and under the 
same hypothesis of Proposition 5.1, we have: 

„ ., ., E 

where 

7 : = 7 (T) and E := de dT ', 

with d > 1 is a positive constant depending on the initial conditions but independent of 
T, and will be given at the end of the proof. 

Remark 6.2 (The constant E depending on time). Let us stress on the fact that, 
throughout the proof and in the rest of the paper, the term E = de dT of Lemma 6.1 might 
vary from line to line. In other words, the term d in the expression of E might certainly 
vary from line to line, but always satisfying the fact of just being dependent on the initial 
data of the problem. The different E 's appearing in different estimates can be made the 
same by simply taking the maximum between them. Therefore they will all be denoted by 
the same letter E. 

Proof of Lemma 6.1. Define the functions u and v by: 

u(x,t) = pt(x,t) and v(x,t) = Kt(x,t). 



= 70/2, 70 e (0,1). 
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We write down the equations satisfied by u and v respectively: 

' ut = (1 + e)u xx - tv x on V, 

< u\ St = 0, 

, u\ t=0 = u° := (1 + e)p° xx - tk° x on /, 

■ Pxx . j -) j~)Pxx 

vt = ev xx H u x + Bu xx - B v x - tu x 



(6.2) 



v\s T = °> 

v\ t =o = f u := en xx 4 



_ ^,_0 , PxPxx _0 



on V, 



on /. 



(6.3) 



The proof could be divided into three steps. As a first step, we will estimate the L°°{V) 

2 1 

norm of the term v x = Kt x . In the second step, we will control the W 2 ' (T>) norm of 
v = n t . Finally, in the third step, we will show how to deduce a similar control on the 
W^' 1 ^) norm of k xx . 

Step 1. (Estimating H^Hoo/d) 

It is worth recalling the equation satisfied by k: 

p x ( ) ' < 

Kf = £K XX H 



- rp x 



In Proposition 4.3, we have shown that k G C 3+a ' 2 . Therefore, writing the parabolic 
Holder estimate (see (2.4)), we obtain: 



I 11 1 |(3+a) 

| K tx||oo,'D — | K lx> — c 



H 





PxPxx 




Ux 



(1+a) 



+ 1/9 

T \Px\d 



(6.4) 



where the term 1 comes from the boundary conditions, and > is a positive constant 
that can be estimated as c H < E (see Remark 2.2). We use the elementary identity 



\ ( v +a) 
to the term 



px pxx 
K>x 



I _u „|r|U+") 

with / = and g = p xx , we get: 



(1+a) 
V 



I (a) 



73 



< 3| P |g +a) 



\Pxx 1 1 00,13 \ / 4 \\Pxxx ||oo,13 \ 



K : , 



x I V 



f'':I 



X I V 



2\p\ (2+a) 

Pxx 00.13 4 Kra oo.p)) 

7 



(6.5) 



where we have used the fact that k x > 7 and k x > \p x \. We plug (6.5) in (6.4), we 
obtain: 



^< E [i + \p\^ + (^Y +a) + \p\^(^ (,,: 



x I v 



x I x> 



4 - ( 1 - |/v 

7 



(2+a) 
V 



(6.6) 
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< E. 



where we have used used the fact that the term has an exponential bound (see 

Proposition 5.2) of the form 

/ / \ (l+a) 

Step 1.1. f Estimating {^J 



From the definition of the Holder norm, we see that in order to control \%r-)_ , it 
suffices to control the three quantities: 



(l+a) 



v 



Pxy * ' liPx 

, K x/t,V \\ K x/ x l X ,V 

We use the the following identity: 



(a) 



and 



Px 



x/ xl t,V 



(§) 



(a) 



LX < 

9/ t ,v 



oo,2? 



oo,T> 



isti + 



oo,T> 



with f = p x and g = n x , we get 



- ~ ( (Px)t,V + ( K x)t,V ' 



xl t,V 



(6.7) 



Similarly, we obtain: 



Pa; a; 
K. 



(a) ||Pxx||oo,Z> / ,(a) 
- -,2 \ K x)x,V + 



x/D 



s / x,X> 



7" 



7 



(6.8) 



We also use the inequality: 



(fg) x %<\\f\\ooA9)xt + \\9\\ooMf)x% 



with / = ^f* and o = ^, we get: 



\ ( Q ) / V n ■ n : 

l^xxPx \ ^ \ xx lx,T> \\^xx\\oo,T> I \(p) || ^xx || oo,2? / 



< 



v x / 1,25 



7 



+ 



-<P*>x% + 



(6.9) 



Similarly, we get 



and 



f«\ , .(f) 

Pxx ||pxx|| o,P / v (f) , (Pxx) t ,p 

^ " \ K x)t,V + 



ft 



a: / t,x> 



(6.10) 



HxxPx \ ^ 2 ^ , ( K *x) 



< 



(f) 
t,2? 



v x / t,£> 



7 



|| K xx ||oo,X> / \(f) _|_ ll K a;x ||oo,X> ; \(f) 



(6.11) 



Collecting the above inequalities (6.7), (6.8), (6.9), (6.10), and (6.11) yield: 



P*\ (l+a \ E fi_u |(2+«)^ N „ 

— > < ~2 ( 1 + + ll«**l|oo,l> 

Kx I j) 7 v 



(6.12) 
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where we have used the fact that 1 < ^, 7 < 1 and \p\^, +a ^ < E (see Proposition 5.2). 



Step 1.2. ( Estimating |/>|g +a) and |ft|g +a) ) 



Using Holder estimate for parabolic equations (estimate (2.4) of Proposition 2.1), and 
similar computations to that of the previous step, we deduce that: 

E 



i(2+a) 



<f(i + w$?), 



and 



Step 1.3. (The estimate for Halloo,©) 



(a) 



By combining (6.6), (6.12), (6.13), (6.14), and by using the fact that \k x \d 
exponential estimate (see estimate (5.1) of Proposition 5.1), we deduce that: 



■ .(3+a) < E_ 

which will be useful later, and as a particular subcase, we have: 

E 

W^tx 1 1 OO, T> ^ ^3 ! 

where we have frequently used that 7 < 1, and we have always taken the maximum of 
all the exponential bounds of the E = de dT form. 



(6.13) 
(6.14) 

has an 
(6.15) 

(6.16) 



Step 2. (Estimating ||Kt|| w 2,i^) 
Step 2.1. (Estimating |M| w 2,i^) 

We use the L 2 estimates for parabolic equations (Theorem 2.3) to the function u satisfying 
(6.2), we obtain: 

IMI^i^ < E(l + Kl|2,x>)- (6.17) 

The term 1 in (6.17) comes from estimating the initial data u°. Since v x = K tx , we plug 
the estimate (6.16) obtained in Step 1.3 into (6.17), we get 



E 



(6.18) 



Step 2.2. (Estimating 11^11^2,1^) 



Arguing in a similar manner as in the previous step, we obtain the following estimate 
for the function v, the solution of the parabolic equation (6.3): 



\v\\ w 2,i (v) < E ( 1 + 



p.r.r 



+ \\B\ 



co,T> 



oo,V 



\Ux\\2,V + ||-B|| 00 ,X)||«xx||2,X> 



oo,T> 



l^xll 2,7? + INxlb,© , 



(6.19) 
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and hence, from (6.16), (6.18), and doing some computations, we deduce from (6.19) 
that: 

E 

The goal of Step 2 follows since v = Kt- 
Step 3. (Estimating ||re xx 11^,2,1^) 

The estimate of | |re x:E | 1^2,1^ requires a special attention. We will mainly use the equa- 
tions on p and re. The four parts || || 2,x>, ||«xxt||2,x», ||«xxx||2,x> and ||« xxxx ||2,x> of the 
above norm will be estimated separately. 

Step 3.1. (Estimate of ||k X x||2,x>) 

Inequality (6.13) directly implies that 



E 

oo,X> — j 

7 



(6.21) 



hence |k xx || 2 ,x> < f . 



Step 3.2. (Estimate of ||re xxxx ||2,x>) 



We first derive the equation on p two times in x, we deduce (using (6.18)) that ||/9 xxxx ||2,x> 
has the same upper bound as ||k xxx ||2,d, i.e. 



E 

\Pxxxx\\2^D ^ ^3- 



(6.22) 



We derive the equation on re two times with respect to the variable x, we obtain: 

, ^PxxPxxx ^xxPxx , PxPxxxx PxPxxx^xx 

^txx — £K>xxxx H 9 r q 

tc IC IC IC 

rv x i\ x rv x rv x 

Pxx^xx PxPxx^xxx Px^xxPxxx ^^xxPxPxx 

9 n 9 I 3 TPxxxi 

tC IC IC tC 

and we use (6.22) and our controls obtained in the previous steps, in order to deduce 
that: 

E 

K^xxxx 2,D — 7- 

T 

In fact, the highest power comes from estimating the following term: 

< 



2 

K xxPxPxx 



2,X> 



K xxPxx 



oo,V 



T 



where we have used the L°° estimate of ||re X x||oo,x»- All other estimates are easily deduced. 
Let us just state how to estimate the other term were ||re X x||oo,£> interferes. In fact, we 
have: 



PxPxxx^xx 



h 7 



< 



2,U 



oo,T> 



\Pxxx\\2,V < — • 



E 

ry3 
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Step 3.3. (Estimate of ||re xx t||2,x> and \\k xxx \\ 2 ,v) 
As an immediate consequence of (6.20), we get 

E 

\\Kxxt \\2,V < —a- 

7 

Deriving the equation on k with respect to x, we obtain: 

2 

, Pxx , PxPxxx Px K xxPxx ,„ 
Ktx — EK-xxx ~r T 9 T~Pxx- {O.Z,o) 

The estimate (6.14) gives 

E 

\\Pxxx \\2,V < — j 

7 

which, together with (6.21) and (6.23), give ||k X ex||2,z> < We deduce as a conclusion 
that: 

, E 

\\ K xx\\ w ^fp) < y , 

and this terminates the proof. □ 

We move now to the main result of this section. 
Lemma 6.3 (W^ 1 bound for p xxx ) 



Under the same hypothesis of Lemma 6.1, we have: 



- E 

\\Pxxx\\ w ^(p^ S 

Proof. Set 



Pxxx\\ W ^(p) < y ■ (6.24) 



r r 
(1 + e)^ 1 + e 

and 

W = Pxxx ~ K - 

We write down, after doing some computations, the equation satisfied by w: 

T 

w t = {1+ e)w xx - ^ +£ y Ktt on V 

w x \ St =0 on St (6.25) 
,„l ._ 7 „o_ n o + 2g ) n T P°xP°xx | ^ 2 n 



(1 + e) 2 ^ (1 + e) 2 Kg (1 + e) 2 ^' 

Here w x \s T = can be checked by deriving the equation satisfied by p with respect to x 
and then with respect to t, and by using the equality (1 + e)p xx = tk x satisfied on the 
boundary dl (which is a consequence of the compatibility conditions). Applying the I? 
theory with Neumann conditions (see for instance [21, Chapter 4, Section 10]) to (6.25), 
we get that p xxx = w + R satisfies 

\\Pxxx\\ w 2,l(p) < E (l + || /e tt || 2 ,z> + \\K t \\ w 2,l {v) + ||«xx ll^.lp)) , (6.26) 
and eventually (6.26) with Lemma 6.1 gives immediately the result. □ 
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7 An upper bound for the BMO norm of p xxx 

This section is devoted to give a suitable upper bound for the BMO norm of p xxx . This 
result will be a consequence of the control of the BMO norm of a suitable extension 
of k xx . The goal is to use this upper bound in the Kozono-Taniuchi inequality (see 
inequality (2.13) of Theorem 2.13) in order to control the L°° norm of p xxx . We first 
give some useful definitions. 

Definition 7.1 (The "symmetric and periodic" extension of a function). Let 

f G C(It), we define f sym (constructed out of f) ouerRx (0,T), first by the symmetry of 
f with respect to the line x = over the interval (—1,0), and then by spatial periodicity. 

Definition 7.2 (The "antisymmetric and periodic" extension of a function). 

We define the function f as v m in a similar manner as f s v m , where we take the antisym- 
metry of f instead of the symmetry. 

We start with the following lemma that reflects a useful relation between the BMO norm 
of f sym and f asym . 

Lemma 7.3 (A relation between f s y m and f as v m ). Let f G C{Lf), then: 

\\f Sym \\BMO(Rx(0,T)) < c (ll/ aSJ/m ||BAfO(Rx(0,T)) + m 2Ix(0,T) (l/^™!)) > 

where c > is a universal constant. 

The proof of this lemma will be presented in Appendix B. The next lemma gives a control 
of the BMO norm of (n xx ) asym . 

Lemma 7.4 (BMO bound for (n xx ) asym ). Under hypothesis (HI), and under the 
same hypothesis of Proposition 5.1, we have: 

\\{n xx ) asym \\ B MO{^T))<ce cT , (7.1) 

where c > is a constant depending on the initial conditions (but independent ofT). 
The function {n xx ) asym is given via Definition 7.2. 

Proof. Let R,(x,t) = n{x,t) — k°(x). We notice that k\s t = 0, therefore K as v m satisfies: 

( \asym asym 

R a t sym - sRZ ym = KPx [ Kx) JyZ ~ T(Px) asym + e(K xx T sym on Rx(0,T) 
R asym {x,0) =0. 

We already know that the right hand side of (7.2) is bounded in L°° by E = ce cT , and 
hence (using Theorem 2.10) the result follows. □ 

We now present the principal result of this section. 
Lemma 7.5 (BMO bound for 

Pxxx)- Under the same hypothesis of Lemma 7.4, we 

have: 

\\Pxxx\\bmo{V) < E. (7.3) 
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Proof. Take 

TK 

1 1 + £ ' 

the equation satisfied by v reads: 

/ -i | \ £T T PxPxx T 2 

1 + £ 1 + 1 + £ 

vlt=o = ^° := pi- — K ° on I 

It ^ 1 + £ 

where we have used the compatibility conditions to check that v x \s T = 0. We can assume, 
without loss of generality, that the initial condition v = 0. This is because being non- 
zero just adds a constant depending on the initial conditions in the final estimate that 
we are looking for. From the fact that v x \s T = 0, we can easily deduce that the function 
v sym sa ti s fi es: 



t 2 r (n ) s y m (n ) s y m ft 
sym _ /, , \ sym , J_ / \sym _ 1 KPx) \Pxx) J—( K Y'^ 

on Rx (0, T) 



^v sym (x,0)=0 on R, 

therefore, using the BMO estimate (2.11) for parabolic equations, to the function v, one 
gets: 

ll v S m ||BMO(Rx(0,T)) < c[\\g\\ B MO(Rx(0,T)) + m 2Ix (0,T) (M)] , 

where c = c(T\) > with < T\ < T. From Propositions 5.1, 5.2, we deduce that 

\\9\\bMO(Rx(0,T)) <E+ \\{Kxx) Sym \\BMO(Rx(0,T)), 

and 

m 2Ix(0,T)(\g\) < E + m2 Ix ( 0>T )(\(Kxx) Sym \)- 

Recall the definition of the term E from Remark 6.2. At this stage, we write the following 
estimate: 

||(^ :E ) SJ/m || J BMO(Rx(0,r)) < C [\\(K X x) aSym \\BMO(Rx(0,T)) + ™2Jx (0,T) (I {^xx) Sym \)] , (7.4) 

which can be deduced using Lemma 7.3. The constant c > appearing in (7.4) is 
independent of T. Finally, we deduce that: 

\\ v x V x m \\BMO{Rx{<d,T)) < c (E + T~ 1 / p \\K xx \\ Pt T}) ■ 

From (6.1), (5.11), we know that r _1 / p ||K a;a ;||p i © < T^ 1 ^ P E. From the previous two 
inequalities, and since v xx = p xxx — ^jff, we easily arrive to our result. □ 
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8 L°° bound for p. 



XXX 

In this section, we use the results of Sections 5, 6 and 7, in order to give an L°° bound 
for p xxx via the Kozono-Taniuchi inequality. 

Proposition 8.1 (L°° bound for 

Pxxx)- Under hypothesis (HI), and under the same 
hypothesis of Proposition 5.1, we have VT > 0: 

\\Pxxx\\ooJ T <^(l + log+^ . (8.1) 

Proof. For T <T\ = T*, where T* is the short time existence result (see Theorem 4.2) 
given by (4.4), inequality (8.1) directly follows. In the other case where T > T\, we apply 
the parabolic Kozono-Taniuchi estimate (2.13) to the function p xxx , together with (7.3) 
and (6.24). Remark that the || Pxxx Hl^/t) can ^ e easu y estimated by the term E. □ 

Proposition 8.2 (A priori estimates) . Under the same hypothesis of Proposition 
3.1, the solution (p,k) G C 3+a '^ (It) satisfies for every < t < T: 

e b(t+l) 

K x (.,t)>e- e >0, (8.2) 

\p(;t)\r> < e e and \k(., t)\f +a > < e e . (8.3) 

Here b > is a positive constant depending on the initial conditions and the fixed terms 
of the problem, but independent of time. 

Proof. Remark that if we consider a function 7 satisfying (3.2), then the right hand 
side of (3.2) can be estimated using (8.1) as follows: 

-(co + ||pxxx(-,OIU-(/))>-^(l + |log7(i)l), E = E{T) = de dT . (8.4) 

This is the motivation to consider the solution 7^ of the following ordinary differential 
equation: 

f7; = -£(l + |log 7T |)7 T , t€(0,T) 
\7t(0)=7oA 

where 70 is given by (3.1). Then, using a continuity argument, joint to the fact that 
fn(t) > 7 2 (t) (see Proposition 3.1), it is easy to check that both (8.4) and (3.2) are 
satisfied with 7 = 7t- Let us now define 

7(T) := 7T (T) < 7r (i) for t G [0, T]. 

Then we deduce from (3.3) that 



k x (.,T) > V9(T) + { Px {.,T)Y > 7(T), VT>0, 

and finally, solving (8.5) explicitly, inequality (8.2) directly follows. Therefore, from 
(6.14) and (6.15), we easily deduce (8.3). □ 
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9 Long time existence and uniqueness 



Now we are ready to show the main result of this paper, namely Theorem 1.1. 

Proof of Theorem 1.1. Define the set B by: 

B - { T > °5 31 solution (j>, «) G C 3+a '^(7r) of 1 
~ \ (1.1), (1.2) and (1.3), satisfying (1.10) J ' 

This set is non empty by the short time existence result (Theorem 4.2). Set 

Too = sup ,6. 

We claim that = oo. Assume, by contradiction that < oo. In this case, let 
5 > be an arbitrary small positive constant, and apply the short time existence result 
(Theorem 4.2) with T = — 6. Indeed, by the tri-exponential bounds (8.2) and 
(8.3), we deduce that the time of existence T* given by (4.4) is in fact independent of S. 
Hence, choosing 5 small enough, we obtain To + T* G B with Tq + T* > and hence a 
contradiction. □ 



10 Appendix A: miscellaneous parabolic estimates 

Al. Proof of Lemma 2.4 {LP estimate for parabolic equations) 

As a first step, we will prove the result in the case where e = 1, and in a second step, we 
will move to the case e > 0. It is worth noticing that the term c may take several values 
only depending on p. 

Step 1. (The estimate: case e = 1) 

Suppose e = 1. Since u = on dl x [0, T], we take u = u asym (see Definition 7.2). Also 
consider the function / = f as v m . Define u by 

u = 

with 

(f) n (x) = 1 if x € (0,2n) 
n ( x ) = O if x > 2n + 1 or x < -1. 

This function satisfies 

u t = u xx + f, on Rx(0,T) 
u(x,0)=0, on R, 

with 
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The proof that 

H^tllp^T \\ u xx \\p,It < 4f\\ p ,i T (10.1) 
can be easily deduced by applying the C alder on- Zygmund estimates to the function u 

2 1 

satisfying the above equation, and passing to the limit n —> oo. Now, since u G W p ' (It) 
with u\t=o = 0, we use [21, Lemma 4.5, page 305] to get 

\\u\\p,i t ^ cT{\\ u t\\ P ,i T + \\u xx \\ p ,i T ) (10.2) 

and 

||wx||p,/ T ^ cVT(\\u t \\p t i T + ||«xx||p,/ T )- (10.3) 
Combining (10.1), (10.2) and (10.3), we deduce that 

jdMlp^T + \\ U x WpJt + ll u xx||p,/ T + ll n *llp,/T — C ll/llp,-fT- 

Step 2. (The estimate: general case e > 0) 

To get the general inequality, we consider the following rescaling of the function u: 

u(x,t) = u(x,t/e), (x,t) G I £ t, 
which allows to get the desired result. □ 

A2. Proof of Lemma 2.7 (L°° control of the spatial derivative) 

Since u G Wp' 1 (It) for p > 3, we know from Lemma 2.5 that u x G C a,a ^ 2 (lT) for 
a = 1 — |. In this case, we use the estimate (2.8) with 5 = VT, we obtain 

IK||oo,J T < c (p){ t ^ (,\\ u t\\p,i T + \\uxx\\p,i T ) +T%- 1 \\u\\ p j T }. (10.4) 

Remark that the fact that u = 0on the parabolic boundary 8 p It, and that it obviously 
satisfies the equation: 

( u t = u xx + /, with / = u t - u xx 
\u = on d p I T , 

then we can apply estimate (2.7) to bound the term ||it|| P) j T . Hence (10.4) becomes (with 
a different constant c(p)): 

< c(p){T2\\ Ut - u xx \\ p j T +Ti 1 T\\u t - u xx \\ p j T } 

< c(p)T*\\u\\ w 2,i ilT) 

< c{p)T V -W\\u\\ w 2,, (ItV 

and the result follows. □ 
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11 Appendix B: parabolic BMO theory 



Bl. Proof of Theorem 2.10 (A BMO estimate in the periodic case) 

Let / be a bounded function denned on R x (0, T) satisfying f(x + 2,t) = f(x,t). We 
extend the function / to R x R + , first by symmetry with respect to the line {t = T} and 
after that by time periodicity of period IT. Call this function /. Set u as the solution 
of the following equation: 

I u t = eu xx + f on R x R+ 

I xx J + 11 ^ 

[u{x,0)=0. 

We apply the standard result of BMO theory for parabolic equations. Since 
/ G L°°(R x (0,T)), then / G BMO(R x R + ), and hence we obtain 
that Ut,u xx G BMO(M x R+), with the following estimate: 

H^llsMO(RxE + ) + ll^xx 1 1 BMO (RxR+) ^ c \\ f \\ B M O (SlxR + ) i 

and hence (from the definition of the BMO space) , 

H^t||sAfO(Rx(0,T) + \\Uxx\\bMO(Rx(0,T)) < c II/IIbMO(RxR + )- (H-3) 

The BMO theory for parabolic equations, particularly estimate (11.2) is rather classical. 
This is due to the fact that the solution of (11.1) can be expressed in terms of the heat 
kernel T defined by: 



,2 

I , i n-i/^. — 

T(x,t) 



(47ret)" 1/2 e"tt, for t>0 
for t < 0, 



in the following way: 

u(x,t)= [ T(x-Z,t-8)f(Z,s)dSds. 

JRxR+ 

As a matter of fact, it is shown in [10] that T xx is a parabolic Calderon-Zygmund kernel 
(here we are working in nonhomogeneous metric spaces in which the variable t accounts 
for twice the variable x). Therefore T xx : BMO — > BMO is a bounded linear operator. 
This result is quite technical and can be adapted from its elliptic version (see [2, Theorem 
3.4]). It is less difficult to show that T xx : L°° — > BMO, a bounded linear operator (see 
for instance [15, Lemma 3.3]). 

Having (11.3) in hands, it remains to show that 

II/IIbMO(RxE+) < c {\\f\\BMO(Rx(0,T)) + m 2/x(0,T) 

with c > independent of T. This can be divided into two steps: 
Step 1. (Treatment of small parabolic cubes) 
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We consider parabolic cubes Q r = Q r (xo,to), {xo,to) € R x R + , with r < y/T. Let us 
estimate the term 

777 / \f~ m Qrf\- 

Assume, without loss of generality, that T < to < 2T. In fact, any other case can be 
done in a similar way because of the time symmetry of the function /. Two cases can be 
considered. If r 2 < to — T then the cube Q r lies in the strip R x (T, 2T) and in this case 



77 / \f-m Q J\ < ||/||bmo(Mx(o,t))- 



IQr 

The other case is when r 2 > to — T. In this case, define and Q£, the above and the 
below parabolic cubes, as follows: 

Qr = Q r (x ,T + r 2 ) and Q b r = Q r (x ,T). 

Since 

T - r 2 < t - r 2 < T < t < T + r 2 , 
then Q r C (Q£ U Q b r ). Moreover, we have \Q r \ = \Q%\ = \Q b r \. We compute: 

mIj f - mQjl - mL lf - 2mQ ' i+mQ - fl 



We remark (from the symmetry-in-time of the function /) that rriQaf = rrigbf, and 

/ \f~c\ = [ |/-c|, VcGR. 

JQr J Qr 

Therefore the above inequalities give: 

777 / \f~ m Q r f\ < 16||/||bmo(rx(o,t))- 
Step 2. (Treatment of big parabolic cubes) 



Consider now parabolic cubes Q r C R x R + , r > \JT. Suppose first that r > 1. Because 
of the symmetry-in-time of the function /, and its spatial periodicity, we compute: 

1 f ~ ~ 2 f ~ 2N f 

777 / \f- m Qrf\<T7T\ l/I^Tol/ I/I' 

IVrl JQ r IVr| 7Q r |Vr| J2Ix(0,T) 

where A is the minimum number of domains D of the form D = (k, k + 2) x (nT, (n+l)T), 
G Z and n G A, that cover Q r . Here 

\Q r \ ~ A x |2J x (0,T)|, A>1. 
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Therefore, the above inequalities give: 

7777 / \f ~ m Q r f\ <cm 2Ix ( ,T)(\f\)- 

IVr| JQ r 

Now suppose that VT < r < 1. In this case we use the fact that < Xi < T, we 
compute: 

,757/ \f-m Q J\<^-[ \f\<^[ \f\<^Lf 1^1- 

IVr-l jQ r |Vr| jQ r |^r| J2Ix(Q,T) T{' J2Ix(0,T) 

Here < and hence 

1 r 

\f-mQrf \ <c( T i)m 2 i x (Q,T){\f\)- 



\Qr 

Steps 1 and 2 give the required result. □ 
B2. Proof of Lemma 7.3. We divide the proof into two steps. 
Step 1. (Treatment of small parabolic cubes) 

Let us consider parabolic cubes Q = Q r (xo,to) C R x (0, T) with < r < ^. Assume, 
without loss of generality, that 1 < xq < 2 (the other cases can be treated similarly). 
Define the left and the right neighbor cubes of Q r {xo,to) by Q~ = Q^{^ — r,to), and 
Q + = Qt(l + r '*o) respectively. Since 2r < 1, then 

Q~ C (0, 1) x (0, T) and Q+ C (1, 2) x (0, T). 

Using the fact that for any function g 6 L x (0): 

/ |s-mn(s)| <2 / | 5 -c|, Vc € R, 

We compute: 

^J Q \f sym -m Q (m\ < ^j Q \f sym +™ Q +{r ym )\ 

< j^rJ Q _\f sym + m Q+ (f asym )\ 

+ i^f Q+ \f sym + rn Q+ (r ym )\- (11-4) 

We know that from the properties of f sym and p sym that m Q+ (/ a ^ m ) = -m Q -{f sym ), 
and 

y S J/m = _ f asym Qn q+ &nd fS ym = jasym Qn q- 

Using the above two inequalities in (11.4), we get: 

W\ S Q lfSVm ' mQ{fSVm)l - ^f^BMOiR^T))- 
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Step 2. (Treatment of big parabolic cubes) 

Consider parabolic cubes Q = Q r C K x (0, T) such that r > \. In this case, we compute: 

1^1 / \f sym - m Q {r m )\ < A / \r m \ < ^ / im, 

M JQ M JQ \Q\ J2Ix(0,T) 

with 

\Q\~Nx |2/ x (0,T)|, 

therefore 

l^i / ir ym - ™ Q (/ s m < cm 2/x(0)T) (i/^D, 

where c is a universal constant. Steps 1 and 2 directly implies the result. □ 
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